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Oh We analyze the statistics of an estimator, denoted by £t and referred to as the slave, for the 

equilibrium susceptibility of a one dimensional Langevin process xt in a potential 4>{x) . The 
susceptibility can be measured by evolving the slave equation in conjunction with the original 
Langevin process. This procedure yields a direct estimate of the susceptibility and avoids the 
' need, when performing numerical simulations, to include applied external fields explicitly. The 

success of the method however depends on the statistical properties of the slave estimator. The 
joint probability density function for x t and £ t is analyzed. In the case where the potential of the 
' system has a concave component the probability density function of the slave acquires a power law 

tail characterized by a temperature dependent exponent. Thus we show that while the average 
value of the slave, in the equilibrium state, is always finite and given by the fluctuation dissipation 
relation, higher moments and indeed the variance may show divergences. The behavior of the 
power law exponent is analyzed in a general context and it is calculated explicitly in some specific 
examples. Our results are confirmed by numerical simulations and we discuss possible measurement 

■ discrepancies in the fluctuation dissipation relation which could arise due to this behavior. 

e 

. I. INTRODUCTION 

a 

o 

U ' A standard experimental technique for probing a system is to measure its response to a small external field. In 
equilibrium, static response functions are related through the fluctuation dissipation relation to appropriate static 
correlation functions. A way to measure such responses in the context of numerical simulations is the slave equation 
method which is used for Langevin type systems, and more precisely in the context of stochastic quantization Via 
qq i the static fluctuation dissipation theorem the slave equation method can be used to compute correlation functions and 
(T) ' has been successfully exploited in numerical simulations of quantum field theories and statistical spin systems 0, . 
I/"") , The advantage of the slave equation method is that it provides a way of measuring cumulant correlators directly with 
properly estimated statistical errors. As has been noted however @, [|, there are circumstances in which the slave 
equation method breaks down. These difficulties arc clearly illustrated by the simple model investigated in this paper. 
Another advantage of the slave method is that the response can be measured without imposing a small external field 
and thus unperturbed correlation functions can be measured simultaneously. Recently various numerical methods for 
measuring response functions in discrete spin systems, without applying an external field, have been proposed . 

Although not the subject of this paper, it is worth noting that for systems in equilibrium more general fluctuation 
dissipation theorems exist that relate dynamical response functions and dynamical correlation functions. In out of 
equilibrium systems that exhibit aging various types of fluctuation dissipation theorems have been to shown to hold 
in mean field models and numerical and experimental evidence points to their validity in finite dimensional systems 
Forms of fluctuation dissipation relations/theorems are also expected to hold in the steady state of certain 
non-equilibrium driven systems. 

The model we investigate is the simplest possible, namely an over-damped particle moving in one dimension in a 
background potential 4>{x) subject to an external field h and thermal noise. The Langevin equation for this system is 



o 
o 



X 



S ; ** = ~(t>'{x t ) + h + r]t , (1) 

where r\t is zero mean Gaussian white noise with correlation function 

(r h r lt ,) = 2TS(t-t') , (2) 

with T the temperature. At zero external field we define the response function £ by 

dxt . , , 

6 = -^|*=o (3) 
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Differentiating Eq. with respect to h and setting h = we obtain the equation of motion of £t as 

& = -4t'{x t )Zt + 1 ■ (4) 

The variable £j is referred to as the slave and Eq.Q the slave equation corresponding to the process Xt . This 
terminology is obvious upon examining Eqs.JQ) and Q) as we see that the process £t is driven by the process xt but 
the evolution of Xt is completely independent of £t . Of course the analysis can be generalized to higher dimensions. 
This is separately interesting and will be addressed in future work. 

Processes whose evolution is similar to that given by Eq. Q arise in a variety of physical contexts such as the 
development of curvature in material line and surface elements and magnetic fields transported by random flows 
H>Q>H>I3' Indeed the behavior exhibited by the probability density functions of the relevant slave variables is exactly 
as analyzed in our simple model. A related slave variable Ct , 

Ct = -5— , (5 

ox 

where xq is the initial value of x t , satisfies the slave equation 

Ct = -<P"(x t )(t . (6) 

Eq. © is identical to Eq. (@J except for the inhomogeneous term +1 on the right. Although the two equations are close 
in form the presence of this inhomogeneous term radically affects the statistics of £t . The behavior of the slave variable 
Ct is is a measure of the sensitivity of a system to its initial conditions and thus related to Lyapounov exponents. It 
was recently shown how Lyapounov exponents could be analyzed via supersymmetric quantum mechanics [Toj . our 
approach is not explicitly developed in terms of supersymmetry but it is clear that it could be rewritten in these 
terms. 

In equilibrium in the presence of an external field h , the statistics of the process x t is described by the Gibbs- 
Boltzmann distribution 

P G B(x) = -^-ex V (-0<i>(x)+phx) , (7) 
Z(h) 

where 

Z(h) = I dx exp(-/3(p(x) + fihx) (8) 



is the canonical partition function for the process x . Clearly one has, by linearity of an expectation of a probability 
distribution, that 

(Oe = -^(x) e , (9) 

where the subscript E denotes expectations taken in the equilibrium state. Using the form of the Gibbs-Boltzmann 
distribution then yields the static fluctuation-dissipation theorem 

(Oe=P((x 2 )e-(x) 2 e ) . (10) 

In simulations, the mean of £ in equilibrium is easy to determine. It therefore provides, up to a multiplicative 
factor, a direct estimator for the variance of x . The numerical usefulness of this result is two-fold. First a straight 
calculation of the variance computed as in Ea. (|10|l is not ideal from the the point of view of precision, as it involves 
the subtraction of two numbers, each potentially much larger than their difference. Measuring £ avoids this difficulty. 
Second, from the variance of £ we have a proper statistical estimate for the error in the susceptibility and the variance 
of x . The utility of the slave method hinges precisely on the existence and size of the variance of £ , the slave variable. 

Given the importance of £ as an estimator to the variance of x by the fluctuation-dissipation relation and its 
direct physical significance as a susceptibility, it is natural to investigate its statistical properties and in particular its 
equilibrium distribution function which we shall denote by p(£) . In this paper, we address this issue in detail. 

The equilibrium probabilty density of £, p(£) , may be obtained from the joint probability density function (PDF) 
P(x, £) via , 



p(0 = J dxP(x,0 



(11) 
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The joint PDF satisfies 13], in equilibrium 



-H FP P( X> + -^ ((cj>" (x)Z-l)P(x,O) = , (12) 
where Hpp is the forward Fokker Planck operator for the process x t and is defined by 

HrrP = - T ^P <V ~ ^'(x)P . (13) 
The Gibbs-Boltzmann distribution for x, P G b(x) of Eq. J7J), is recovered via 

Pgb{x) = [ d£P(x,0 , (14) 



and of course satisfies 

-H FP P GB (x)=0 . (15) 

II. STATIC FLUCTUATION DISSIPATION RELATION 

The fundamental nature of the static fluctuation dissipation relation between the equilibrium susceptibility and the 
systems variance means it is illuminating to verify the relation Eci. (|10f) from Eg. 1121 directly). We set 

P(x,Z)=P GB (x)F(x,0 , (16) 

and define F n (x) , where it exists, by the equation 



F n (x) = J dt£ n F(x,0 . (17) 

Clearly 

F (x) = l , (18) 

and 

(Oe = J dx P GB (x)F 1 (x) . (19) 

It follows from Ea. l(T2)) that 

T^--cl>'(x^-^-F(x,O + -^((<i>' , (x^-l)F(x,O)=0 . (20) 



If we multiply by £ and integrate over all £ we find 



- d>'(x^j ^F x {x) - 4>"(x)Fi(x) + 1 = 0. (21) 



"' ^^-F^x) - cj>'{x)F 1 (x) +x)=0 , (22) 



This can be rewritten in the form 

dx \ dx 
which implies that 

T-^-F 1 (x)-(t) , {x)F 1 (x)+x-a^0 , (23) 
ox 

where a is an integration constant. If we then multiply by Pcb{x) and integrate over all x we find 

a=(x) E . (24) 
We now multiply by Eq. (|23|l by (x — cl)Pqb(x) and integrate over all x to obtain 

-T I ' dx P GB (x)F 1 (x) + { dx P GB (x)(x - a) 2 = , (25) 



which is precisely the static fluctuation dissipation relation Eq. I|10|) . The derivation therefore confirms the properties 
of the estimator as deduced from the stochastic differential equation for Xt ■ 
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III. GENERAL PROPERTIES OF THE ^-PROBABILITY DISTRIBUTION 

It follows from Eq. @ that for £ ~ , £t ~ 1 . This implies that there is always a positive flow of probability from 
negative to positive £ . In equilibrium therefore the support for P{x 1 £) lies in range £ > . We shall assume also, 
which can be justified subsequently, that P(x, 0) = . 

The Laplace transform of the joint probability function is 

/>oo 

P(x,s) = / d^e- si P(x,£,) . (26) 



It follows from our discussion that we expect P{x, s) — > faster than s 1 as s — > oo , and 



l o 

die-^—P{x,i)=sP{x,s) . (27) 



The large £ behavior of P(x,£) is also of importance. As will become clear later this is strongly influenced by the 
behavior of <j>"{x) . Let there be an interval U such that for x G U , <f>"{x) ~ — g where g > . It follows from Eq.|@} 
while x remains in U then £ will grow exponentially. If the time for which x remains in U is r then the excursion 
experienced by £ will be roughly of the form 

e = £oe sr , (28) 

for some £o . When x leaves U , £ will decay rapidly back to small values. Since the the process x t is essentially 
without memory the distribution of r will be exponential. 

p(r) ~ |ttexp(-Mr) . (29) 

The large £ behavior of p{£) is determined by these excursions. We have 

|) . (», 

where a* = fi/g . Hence we can expect a power law in the distribution for large values of £ when there exists a region 
U in which the potential <fi(x) is concave. Of course the argument above does not allow an easy calculation or even 
estimation of the exponent since the time spent in concave regions is affected by the global structure of the potential. 

A potential which exhibits a region of concavity and for which therefore we expect to find a power law behavior 
for p(£) , is 4>(x) = (1 — £ 2 ) 2 /4 . We have simulated the Langevin process x t and its slave by integrating Eqs.Q and 
using a second order stochastic Runga-Kutta method 12]. Shown in Fig. is the time series obtained for £ at 
(3 = 1/T = 1 . The mean value of the process as predicted by the static fluctuation dissipation relation Eg . 111)1 1 is 
shown by the thick horizontal line. A direct measurement of the time series average confirms this result, as it should. 

The most striking feature of the time series however is that £ t spends most of the time below the average value with 
intermittent spikes rising to large values. These spikes are the means by which the time series fills the power law tail 
in the distribution for p(£) as indicated by the intuitive argument explained above. The implication of this result for 
simulations is that it is essential to include the intermittent upward excursions if correct estimates are to be obtained 
for the susceptibility. If therefore one were to measure (£}e by taking a time series average, then in order to obtain a 
satisfactory estimate one must ensure that the time interval of the measurement is of a length sufficient to sample the 
rare but large excursions that represent the power law tail of p(£) . The typical value of £ that one measures during 
a simulation is well below the mean value. Thus measuring over too short a time scale, or for other reasons omitting 
the large excursions will lead to one to an under estimate of (£)e ■ 

Another approach to the fluctuation dissipation theorem would be to be to compute numerically or measure phys- 
ically the expectation value £ a of the slave variable £ and the variance of the the original variable x and use these 
results to provide an estimate T e f / of the temperature of the system. 

1 eff = 7 ■ 

If for the reasons discussed above we under estimate £ a then our estimate for the temperature will be too high, 



Teff > T 



(32) 
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T 


a* (numerics) 


0.66 


2.46(3) 


0.50 


1.86(2) 


0.40 


1.67(1) 


0.33 


1.54(1) 



TABLE I: The exponent a* , estimated by straight line fit to the tail of the log-log plot of the numerically generated histogram, 
as a function of (5 for the potential 4>{x) = (1 — x 2 ) 2 /4 . 



It is interesting to note that this result with T e ff/T > 1 is also seen in ag ing systems where the dynamical fluctuation 
dissipation theorem is violated because the system is not in equilibrium 11]. However here the apparent violation is 
due to an error of measurement in an equilibrium state and hence of somewhat different origin. 

The value of the exponent a* appearing in Eq. I|30[) may be obtained from the numerical simulation by fitting a 
straight line to the large £ region of the log-log plot of the numerically generated histogram of £ . We have carried out 
this procedure for a range of values of the temperature T and the results are shown in table (P) . Most importantly 
we see that the exponent a* depends continuously on temperature. In addition, two features of the behavior of a* 
in these numerical results stand out. First, as T — > oo then it appears a* — > oo (as a* increases an accurate fit of 
the power law tail's exponent becomes difficult due to the lack of statistical weight in the tail). Second, as T — ► the 
numerics is consistent with a* — * 1 . This latter result is particularly significant. At low temperature the PDF for x 
takes the form 

Pgb(x) = ~(6(x-1)+5{x + 1)) . (33) 

This implies that (x)e — and (x 2 }e = 1 . It follows from the static fluctuation dissipation relation that 

1 

T 



(Oe =i I , as T — > . (34) 



For large £ we have 



for some positive C . It follows that 

/OO nOO f < si 

warn - J dee^+s - 7-731) • (36) 

Comparing this result with the static fluctuation dissipation relation we see that at low temperature we should expect 
a* — > 1 and more specifically C/(a* — 1) ~ 1/T . This is consistent with the numerical results shown in table (|lj| . 

For a given value of a* the existence of a power law tail for p(£) means that moments (£™}.e diverge for n > a* . 
At low T, where a* < 2, therefore even the variance of the estimator for the susceptibility has become divergent. At 
this point it has ceased to be a useful estimator and provides no reasonable estimate for a statistical error on the 
susceptibility. The slave equation method therefore becomes ineffective under these circumstances. These issues in 
relation to simulations of quantum field theory and spin models in statistical mechanics have been noted before 0, 01 • 



IV. GENERAL THEORY 



The above observations on the nature of the PDF of £ are much clarified and rendered more robust by an under- 
standing of the general theory of the joint PDF. In order to pursue the analysis it is convenient to make a standard 
transformation that renders Hpp into self adjoint form. We define Q(x,£) so that 

PM = QM eM -tt (x)/2) ■ (37) 

It follows that Q obeys 



- H Q(x, £) + — m"(x) - l)Q(x, 0) = 



(38) 



6 



15 



10 



S(t) 



7i. 




200 



FIG. 1: Time series for slave £t with potential cj>(x) = (1 — x 2 ) 2 /4 at /3 = 1 . Shown by the thick horizontal line is its average 
value 



where the manifestly self adjoint operator Hq is given by 



with 



H = -T— + V( X ) 



We introduce the Laplace transform of Q with respect to the variable £ 

/'OC 

Q(a;,s)= / d£ exp(-s£)Q(x, 
Jo 

and thus corresponding Laplace transform of P is given by 

exp (-04>(x)/2)Q(x,a) 



P(x,s) = 



We note that 



P(x,0) = 



Z? 



exp (-(3<f)(x)) 



which thus gives us the initial conditions for Q at s = to be 

exp(-/ty(a0/2) 



Q(x,0) 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 
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The evolution equation for Q is 

H Q + sQ + s^-Q = . (45) 

OS 

We notice that a particular solution of Eq. I|45l) can be written as a power-series expansion in s 

oo 

QW(V) = E S " + "^W ' ( 46 ) 

n=0 

where the indicial equation determining the allowed values of a is 

#o/ (a) =0 • (47) 

The general solution may then be written as a linear superposition of these particular solutions: 

Q(x,a)=J2 w °Q la Hx,s) • (48) 

a 

We notice that a — is always a solution of Eq. (|47f) with a corresponding given by 

/(")(,) = exp(-/30(x)/2) 

The initial condition Eq. 14411 means that m„ = for a < , otherwise Q would diverge at s — . The solution must 
therefore be of the form 

oo 

Q(x,s)=J2s n fi°\x) + J2^Q ia) (x,s) . (50) 

n=0 a>0 

It is easy to see, on reconstructing p(£) from the above Laplace transform, that it must have the form 

p(0 « 7^ ( 51 ) 

for large £ where a* is the smallest strictly positive solution to the indicial equation with w a * 7^ . If there are no 
solutions to the indicial equation with a > then Q(x, s) must be given by 

00 

Q(x,s) = Y,* n fi 0) (x) (52) 

and is analytic in s and all of the moments of £ will presumably exist; thus /?(£) will not have a power law tail at 
large £ . The physical arguments leading to Eq. (|30|l suggest that Eq. (|51(l should hold when there is a region where 
4>"{x) < . In addition the same physical argument and the static fluctuation dissipation relation also show that we 
should have a* > 1 for finite f3 , otherwise £ will diverge which means that the variance of x on the right hand side 
of Eq. I|10|) diverges, which is clearly not possible for a sufficiently confining potential. We shall now confirm this 
physical picture mathematically. If we define the operator 

Ao = VTl + ^(x) , (53) 

then clearly we may write Hq = AqAq which shows that Hq is positive semi-definite. If / is a solution to the indicial 
equation Eq. H47fl . multiplying by / and integrating over all x gives 

J dx fH Q f + ajdx f<j)"{x) = . (54) 

When <p"(x) > for all x , it follows that a < . As expected, no power law behavior is possible in this case. 
The indicial equation Eq. (|47|l may also be written as 

- rJjLf + (tf^ i^x)f - / + in-)f / = , (55) 
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which is equivalent to 



H a f+ a(1 T a) {<j>'[x)ff = Q , (56) 



where the positive semi-definite H a is given by 

H a = A^ a A a , (57) 

with 



This thus yields 



A a = VTl + ^ {x) . (58) 



dx fH a f + a{l T a) I ,/.,-/-(,/(,■) f = n . (r.oi 



which implies that a(l — a) < , thus if a is positive then we must have a > 1 , again confirming the physical 
reasoning of the Introduction. 

We shall now show that if there exists a region where <j>" (x) < , then there is a solution to the indicial equation 
with a > 0, and from the preceding argument, if such an a exits then a > 1 . 

We consider the following eigenvalue equation, 

H *Po + a<f>"(x)ip = E (a)ip , (60) 

where Eq denotes the ground state energy and ip the corresponding ground state wave function. Consider this 
eigenvalue problem at a — , here we have 

M a =o = exp (-/30(x)/2) (61) 
E Q (0) = . (62) 
First order perturbation theory shows us that 

■^E (a) = J dx </> o y'(x) 

= ^|dx(0'(x)) 2 exp(-^(x)) . (63) 

Thus at a — we have -^-Eq{o) > , thus there is a region of a > where Eq(o) > and there can be no acceptable 
solution a to the indicial equation in that region. 
We shall now use the well known variational formula 



Eo(a) = min^, J dx 



T ( + (0'(x)) 2 + 1(1 - 2a)0"(x) j , ^ 



(64) 



where the minimum is taken over all functions such that J dx ip(x) 2 = 1 . Define g = —mm x {4>"(x)} and consider the 
case where g > and let xo be a point where this minimum is obtained. Without loss of generality we take xo = . 
Now consider the trial wave function 

r(x) = (^Y e M~cx 2 /4) (65) 
for c large and positive. Using Eq. (j64(l we obtain 

£o(a)<^+(^(0'(O)) 2 -f(2a - 1)) + f-U . (66) 

Thus if c 3> 1 and a ^> c then we have E (a) < . Assuming the continuity of E (a) , along with the fact that E(a) 
is positive in a region (0, 1) for some I > , we have shown the existence of a* > such that E(a*) = . 
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Now it remains to be shown that in the case where (j>(x) is concave for some range of x , that the coefficient w a * in 
Eq. 148|) is non-zero. We recall the boundary condition P(x,0) = which means that for large s that Q(x,s) must 
decay more quickly than 1/s at large s . 

The Eq. (|45|l may be written as 

Q (x ,s) = -j^G(^,)sl a QV,, W ) , (67) 
where G is the Green's function obeying 

(H + s)G(x,x';a) = 6(x-x') . (68) 
For fixed x and x' and s>lwe have from Eq. ljfi"5jl 

G(x,x';s)n 5 ^- x '^ , ( 69 ) 



which means that for large s 



and hence 



Q{x,s)^-f{x)^-Q{x,s) (70) 
os 



Q(x,s) n H(x)exp(-s/<f>"(x)) . (71) 

This means that in particular if there is a point x where <j) 1S concave then Q^°'(x, s) diverges there and thus the 
full solution needs to have at least one w a ^ as Q^(x, s) has the same divergent behavior there as s — > oo , the 
coefficients must then be chosen to cancel the divergence. A similar mechanism was identified in simplified discrete 
versions of the the slave equation of the model discussed in this paper Q . 

V. SPECIFIC EXAMPLES 
A. The Simple Harmonic Oscillator 

The simplest example one can consider is the simple harmonic oscillator with 

tfs) = ^r • ( 72 ) 

From the theory of the precedent section we know that there should be no power law behavior in p(£) , and although 
the problem can be explicitly solved it is instructive to work through the mathematics as formulated in Section III. 
The indicial equation in this case is 

dx 2J AT J 2 y JJ y ' 

The left hand side of Eq. Q73JI is the Hamiltonian for a quantum simple harmonic oscillator of mass m = 1/2T and 
frequency uj — X . The energy levels are thus E n — (n + i)A and comparing with the right hand side of Eq. I|73() 
immediately yields that the solutions for a are a = — n and hence, as predicted, are all negative. Clearly the slave 
equation reads 

6 = -A£ + l , (74) 
and so the equilibrium distribution of £ is a delta function at the fixed point £ = 1/A . 

p(0 = *(£-i) • (75) 
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B. The Potential 4> = |as| 

We now consider the potential 



(76) 



Again there can be no power law behavior. The interesting point about this potential is that the full distribution of 
£ can be computed. One can easily solve for P(x, s) to obtain 

P{x,s) = ^ew(-H°)\*\)exp(-W*+\) i ) > (77) 



where 



Then after integrating over x , the Laplace transform can be inverted to yield 



P(0 = rr^TTa exp — - - — . (79) 



7ri(/3£)S V 4 ft 

Notice that the large s behavior is not that predicted by Eq. lf7T|l due to the delta function singularity in <j)"{x) at 
x = , where this occurs it is easy to see that one has a behavior of the form exp(— Ay/s) but the conclusions stay 
the same. Indeed the simplest cases exhibiting a power law distribution in £ are those where 4>"(x) is composed of 
delta functions, as is the case for continuous piecewise linear or quadratic potentials <p ■ 

C. The piece-wise continuous quadratic potential 

Here we consider the potential 

0(x) = ^(\x\-l) 2 -hx . (80) 
In this case we again expect a power law tail in the distribution of £ . The indicial equation in this case is 

~ T §^f + -h + lf- \{l - 2a)) / = 0; x < (81) 

along with the continuity of / at x = and the jump condition. 

- T (g| 0+ - g| _) + (1 - 2a)/(0) = . (82) 
The solution which decays as \x\ — * oo is 

fix) = A+D^ a (yp(x - h - 1)) ; .t>0 

fix) = A_D^ a (yfpih - 1 - x)) ; a;<0 , (83) 

where -D p denotes a parabolic cylinder function of index p 15]. Using the jump condition and continuity at x = , 
along with the identity 

TO -fc p (*) + A>+i(*) = - (84) 

we find that a obeys the equation 

T(a,P,h)=0 (85) 
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FIG. 2: Graph of T(a, (3, h) against a at f3 = 1 and h = . Note there is only one strictly positive solution to F(a, j3,h) — . 
where 

T(a,P,h) = 2^(a-l)D- a (-^(l + h))D- a (-y/j3(l-h)) 

- Di_„(- v ^(l-/i))i3-a(->/^(l + fe))-I>i-a(-V^(l + ^))^-o(-Vl9(l-'»)) ■ ( 86 ) 

Numerically solving Eq. (|85() shows that we have one positive root a* , the root a = while all the others are 
negative. An example is shown in Fig 
At small (3 we can show that 



V/3 , 



(87) 



while the large /3 behavior depends on the value of h . In the case \h\ < 1 the system has two local minima and we 
find that as T — > then a* is given by 



1 



1 



[exp(-/3(l + /0 2 /2) + exp(-/3(l - hf /2)] 



Here we see that in the zero temperature limit a* — ► 1 . In addition the asymptotic form for a* Eq. (|88|l tells us that 
the coefficient C of the power law tail must behave as 



C 



£-exp(-/V2) 



(89) 



in the case when h = . This follows from the static fluctuation dissipation theorem and the fact that the variance 
of x is non-zero as T — > when /i — . Thus although the power law exponent a* decreases, the coefficient of the 
power law component of the PDF is, in this case, tending to zero exponentially quickly. This can be understood 
physically as the excursions into the region where <fi is concave are into regions of high energy, whose Boltzmann 
weight is exponentially suppressed. 

In the case \h\ > 1 there is only one local minimum and we find as T — > that 



1 + 



1*1-1 



(90) 
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T 


a* (numerics) 


a* (predicted) 


2.0 


2.54(1) 


2.560 


1.0 


1.61(4) 


1.612 


0.5 


1.22(1) 


1.191 


0.4 


1.12(1) 


1.121 



TABLE II: Exponent a* evaluated from histogram of £ numerically generated for the potential <f) of Eq. I|8()|l regulated at the 
origin with e = 0.05 compared with analytical prediction of Eq. I85H for various values of T and at h = . 

The predictions obtained by our method may be confirmed by numerical simulations. In the numerical simulations 
one needs (/>' to be continuous and so the cusp in the potential at x — needs to be regularized. We take a small 
interval [e, e] where we set 4>'{x) = (e — l)x/e + h , this choice ensures <p'(x) is continuous and in the limit e^Owe 
recover the potential of Eq. JSUJ). Simulation results are performed with e = 0.05 and compared with the analytical 
result for the potential <f> of Eq. Q80JI. A table comparing the exponents obtained from the numerical simulation 
and from solving Eq. (|85|l for the case ft, = is shown in tabic (HQ, the agreement is excellent and the deviation is 
compatible with it being of order e = 0.05 . Similar agreement is found when h is non-zero. 



D. The 'W Potential 

Here we consider the 'W shaped potential given by 



(x) = a\x — 1| 

= bl-i 



\x\<l 
\x\ > 1 



The second derivative of the potential is thus given by 

(f>"(x) = -2aS(x) + (1 + a)5(x - 1) + (1 + a)S(x + 1) , 



(91) 
(92) 



and so from the general theory of section III we expect a power law behavior where a > and a < — 1 . The indicial 
equation is given by 



<9 2 a 2 



(93) 



along with the jump conditions 



2T— /| + o(l-2a)/(0) = 



T 



Tx fW+ ~Tx ni - 



(_-a)(l + o)/(l) 



(94) 
(95) 



The solution of the indicial equation, which decays as \x\ — > oo is 

f[x) = A + exp(/3o|ir|/2) + A_ exp(-/3a|x|/2) \x\ < 1 
= Sexp(-j9|as|/2) \x\ > 1 

Continuity at x — 1 then gives 

B exp(-/3/2) = A+ exp(/3a/2) + A_ exp(-/3a/2) , 

allowing for the elimination of the variable B . The vector 

A+ 
A- 



(96) 
(97) 

(98) 



is then determined by Mu = , where 



M ~ (exp(/3a/2)a(l + a) eatp(-j9o/2)(a(l + a) - a) ) ' (") 

The possible values of the exponent a are then determined by the existence of a solution such that u^O , that is to 
say detM = , which yields the solution a = or 

a = a* = 1 + —. 1 — r . (100) 

(1 + a) (exp(/3a) - 1) V ' 

We see from Eq. <|100[1 that a* is positive and greater than one in the region where has a concave component as 
predicted by the general theory. In the region a > we find that 

a(l + a)/3 
exp(— /3a) 

~ 1 + V as Z 3 00 ■ ( 101 ) 



(l + o) 



In the region a < — 1 we find 



«* ~ 1 + i — \TT~~ — as /? ^ 
|a||l + a|/3 

w 1 + : 1 : as /3 -> oo . (102) 
|l + o| 

We see in the case a < —1 in the limit T — > that a* > 1 , this must be the case from the static fluctuation 
dissipation relation: here there is only one minimum and we therefore have 

R , 2\ „ Pjdx x 2 exp(-/?H|:r|) 2T 

^ )£ * J&expHgHH) ra ^ ' (103) 

and hence (£)e has no divergence as T — > . 

An interesting point emerges here. If the variance of a; is non-zero in the limit T — > the static fluctuation dissipation 
theorem tells us that the coefficient a* — > 1 in this limit. However, even if the variance of x tends to zero then a* may 
still tend to one in the zero temperature limit. This is seen in the case of the potential <p(x) = (\x\ — l) 2 /2 — hx where 
the minima are only degenerate at h = . The average value of £ stays finite because the prefactor, denoted in this 
paper by C , of the power law component of the PDF of £ is tending to zero sufficiently rapidly. In the cases we have 
examined here it seems that a* — > 1 in the zero temperature limit when there are at least two local minima of the 
potential <p ■ It seems possible therefore that the zero temperature behavior of the exponent a* encodes geometrical 
or topological properties of the potential <j) . 



VI. CONCLUSIONS 



We have analyzed the equilibrium distribution of a slave variable £t which is the estimator of the susceptibility 
of a one dimensional Langevin process x% . Even though the equilibrium statistics of xt are such that all moments 
are finite (for a sufficiently confining potential), the probability density function of the slave £ f can have power law 
tails characterized by a temperature dependent exponent. This power law behavior is present when the potential <p 
has a concave component. The origin of this power law can be understood from simple qualitative arguments. The 
behavior of the power-law exponent can be analyzed in the steady state using the Fokker-Plank equation for the 
equilibrium joint probability density function P(cc,£) for (xt,£,t) ■ The exponent is large at high temperatures and 
decreases on decreasing the temperature. As the temperature is reduced its higher order moments of the slave £t 
diverge and ultimately the variance of the slave may diverge thus rendering it a poor estimator for the susceptibility. 
This pathology in the slave statistics had been observed in Langevin simulations of spin and quantum systems 0, |j| • 
A number of exactly soluble cases were analyzed and the results confirmed by numerical simulation. 

In future work it would be interesting to generalize our results to higher dimensional systems, notably interacting 
systems where phase transitions may occur. The temporal evolution of the PDF of the slave is also worthy of future 
study. It would be interesting to know how quickly the tails of the slave's PDF fill out and after what time it becomes 
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equilibrated. One would also like to understand over which timescale temporal averages need to be carried out in 
order to numerically verify the static fluctuation dissipation relation. Finally the analysis developed here could prove 
useful in the analysis of similar slave variables occurring in Langevin systems. 
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